#!/usr/bin/env Rscript

############################################################
### Least squares estimation
### y = beta0 + beta1*x --> Y = (1, X)beta
### Y = XB ---> X 不一定可逆
### t(X)Y = t(X)XB -->  inverse(t(X)X)t(X)Y = B
############################################################

X = cbind(rep(1, 5), c(1.00, 2.00, 3.00, 4.00, 5.00))
Y = c(3.70, 4.20, 4.90, 5.70, 6.00)

beta.hat = solve(t(X) %*% X) %*% t(X) %*% Y
